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Oxygen isotope effects on the transport properties have been studied in high-quality epitaxial 
thin films of Lao.75Cao.2sMn03 and Ndo.7Sro.3Mn03. In the paramagnetic state, the resistivity 
can be well fitted by p(T) = (A/s/T) exp(E p /kBT) with the parameters A and E a depending 
strongly on the oxygen isotope mass. The resistivity below 80 K almost perfectly follows p — 
p + Blo s / sinh 2 (huj s /2kBT) with Tiui s /kB ~ 100 K. Both p and B increase by about 15(3)% 
upon raplacing le O by ls O. The results provide quantitative constraints on the basic physics of 
manganites. 



The discovery of "colossal" magnetoresistance (CMR) 
in thin films of Rei^A^MnOa (Re = a rare-earth ele- 
ment, and A = a divalent element) Q has stimulated 
extensive studies of magnetic, structural and transport 
properties of these materials [||. The physics of man- 
ganites has primarily been described by the double- 
exchange (DE) model However, Millis, Littlewood 
and Shraiman pointed out that the carrier-spin inter- 
action in the DE model is too weak to lead to the carrier 
localization in the paramagnetic state, and thus a second 
mechanism such as small polaronic effects should be in- 
volved to explain the observed resistivity data in doped 
manganites. The central point of the model is that in the 
paramagnetic state the electron-phonon coupling con- 
stant A is large enough to form small polarons while the 
growing ferromagnetic order increases the bandwidth and 
thus decreases A sufficently to form a large polaron metal- 
lic state. Many recent experiments PHldfl have provided 
strong evidence for the existence of small polarons in the 
paramagnetic state, and qualitatively support these the- 
oretical models 

Alexandrov and Bratkovsky j?J have recently argued 
that the model suggested by Millis et ai, cannot quanti- 
tatively explain CMR since the characteristic theoretical 
field (~ 15 T) for CMR is too high compared with the ex- 
perimental one (~4T) (|| . They thus proposed an alter- 
native theory for CMR. The basic idea of their model is 
that the small polarons form localized bound pairs (bipo- 
larons) in the paramagnetic state due to strong electron- 
phonon coupling while the competing exchange interac- 
tion of polaronic carriers with localized spins drives the 
ferromagnetic transition. The transition is accompanied 
by a giant increase in the number of small polarons which 
are mobile carriers and move coherently at low tempera- 
tures. This model appears to be able to explain the CMR 
quantitatively. 

In order to discriminate between these different mod- 



els, we study the oxygen isotope effect on the trans- 
port properties of high-quality epitaxial thin films of 
Lao.75Cao.25Mn03 and Ndo.7Sro.3Mn03. The intrinsic 
resistivity of these compounds shows a strong dependence 
on the oxygen isotope mass in both paramagnetic and fer- 
romagnetic states. A quantitative data analysis provides 
essential constraints on the basic physics of manganites. 
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FIG. 1. The resistivity of the oxygen-isotope exchanged 
films of (a) Lao.75Cao.25Mn0 3 ; (b) Nd .7Sr .3MnO 3 . 

The epitaxial thin films of Lao.75Cao.2sMn03 (LCMO) 
and Nd . 7 Sr .3MnO3 (NSMO) were grown on <100> 
LaA103 single crystal substrates by pulsed laser depo- 
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sition using a KrF excimer laser Jl5j . The film thickness 
was about 190 nm for NSMO and 150 ran for LCMO. Two 
halves were cut from the same piece of a film for oxygen- 
isotope diffusion. The diffusion for LCMO/NSMO was 
carried out for 10 h at about 940/900 °C and oxygen 
pressure of about 1 bar. The 18 0-isotope gas is enriched 
with 95% 18 0, which can ensure 95% 18 in the 18 
thin films. The resistivity was measured using the van 
der Pauw technique, and the contacts were made by sil- 
ver paste. The measurements were carried out from 5 to 
380 K in a Quantum Design measuring system. 

Fig. 1 shows the resistivity of the oxygen-isotope ex- 
changed films of (a) LCMO; (b) NSMO over a tempera- 
ture range 100-300 K. It is striking that, in all the cases, 
the 18 O samples have lower metal-insulator crossover 
temperatures and much sharper resistivity drop. The 
Curie temperature Tq normally coincides with a temper- 
ature where din p/dT exhibits a maximum. We find that 
the oxygen isotope shift of T c is 15.0(6) K for LCMO, 
and 18.0(6) K for NSMO, in excellent agreement with the 
results for the bulk samples 16 1. 
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FIG. 2. The low-temperature resistivity of the oxy- 
gen-isotope exchanged films of (a) Lao.7sCao.25Mn03; (b) 
Ndo.7Sro.3MnC>3. The solid lines are fitted curves by Eq. 1. 

In Fig. 2 we plot the low-temperature resistivity of 
the oxygen- isotope exchanged films of (a) LCMO; (b) 
NSMO. In both cases, the residual resistivity po for the 
18 O samples is larger than for the 16 O samples by about 
15%. We have repeated the van der Pauw measurements 
at 5 K for several times with different contact configura- 
tions. We checked that the uncertainty of the difference 
in po of the two isotope samples is less than 3%. We 



also found that the isotope effect is reversible upon the 
isotope back-exchange. It is worth noting that the po for 
the 16 O NSMO film is nearly the same as for a single 
crystalline sample jL7|. This indicates that the electrical 
transport properties observed in our films are intrinsic. 

The low temperature resistivity has been explained by 
small polaron metallic conduction [Q, which leads to a 
formula lilJTl: 



p(T) = po + Buj s / sinh 2 (hu s / '2k bT), 



(1) 



where oj s is the frequency of the softest optical mode 
(about 100 K [Q), and B is a constant, being pro- 
portional to m p /n (here m p is the effective mass of po- 
larons, and n is the mobile carrier concentration). The 
residual resistivity po oc m p /riTo [^0|, where K/tq is the 
zero temperature scattering rate which is associated with 
the random potential produced by randomly distributed 
trivalent and divalent cations | pl| . The magnitude of 
%/tq has been found to be 18(3) meV from optical data 
of the compounds Ndo.7Sr .3Mn03 |17| ], Lao.7Sr .3Mn03 
and Lao.7Cao.3Mn03 p2fl . Since the Drude weight (oc 
n/rrip) for Ndrj.7Sr .3MnO3 (l7j is about a factor of 3 
smaller than for Lao.7Cao.3Mn03 p^l , nearly the same 
Ti/tq observed for both compounds implies that the U/tq 
is nearly independent of m p /n, and that both po and B 
are proportional to m p /n. 

The best fit to the data shown in Fig. 2 indeed demon- 
strates that, in both cases, po increases by 15(3)%, and 
B by 18(3). The results indicate that there is a substan- 
tial oxygen isotope effect on m p /n. Since we have shown 
that |l(],[l6| the total hole concentrations for two isotope 
samples should be the same, the observed isotope effect 
on m p /n implies that m p depends strongly on the isotope 
mass. This strongly suggests that the carriers in the fer- 
romagnetic state are indeed small polarons, in agreement 
with the theoretical model Q . We also find that both p 
and B for NSMO are about twice as large as those for 
LCMO, implying that the effective mass of carriers for 
NSMO is larger by a factor of 2.4 than for LCMO. This is 
consistent with specific heat data; the electronic specific 
heat coefficient 7 — 8 mJ/moleK 2 for Lao.sCao^MnOs 
||, while 7-20 mJ/moleK 2 for Nd .67Sr .3 3 MnO3 @. 

In the paramagnetic state, the dominant conduction 
mechanism is thermally activated hopping of adiabatic 
small polarons. The mobility is given by [£5| 

P= e -^^eM-E a /k B T), (2) 

where d is the site to site hopping distance, which is equal 
to a 1 'V2 since the carriers in manganites mainly reside on 
the oxygen sites Q , u is the characteristic frequency of 
the optical phonons, and E a is given by [ p7|j25|j2^ | E a = 
(j)E p /2)f(T) — t, where E p is the polaron binding energy, 
t is the "bare" hopping integral, 77 is a constant (77 = 1 
for Holstein polarons ^5| and rj — 0.2-0.4 for Frohlich po- 
larons §§), and f(T) = [tanh(hu) /4k B T)]/{tiu> /4k B T) 
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If one considers that small polarons can be bound into 
localized pairs (bipolarons) or to the impurity centers 
f20f , one finds [|7],p0| that the mobile carrier density n = 
2(2nm p k B T/h 2 ) 3 ' 2 exp(-A/2k B T) when T<W p /k B and 
for a parabolic band, where W p is the polaron bandwidth 
(irip — 6h 2 /a 2 W p ), A is the bipolaron binding energy |t]] 
or twice the gap between localized impurity levels and the 
bottom of the polaron band ^0|. In fact, the above n(T) 
expression is the same as that for semiconductors when 
the chemical potential is pinned to the mobility edge of 
the bipolarons or to the impurity levels. For a single- 
type (hole or electron) carrier system, A = 2E S , where 
E s is the activation energy deduced from thermopower 
data. Then in the temperature region: Tc<T<W p /kB, 
the resistivity p — 1 /nep is given by 



eME P /k B T), 



(3) 
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and E n 



where A 

E a + A/2. The quantity A should strongly depend on 
the isotope mass because W p = W exp(— TE p /huj ) = 
W exp(— g 2 ). Here W is the bare bandwidth (W = 
12t), r is a constant, which is less than 0.4 for both Hol- 
stein and Frohlich polarons when the coupling constant 
A = 2E p /W o <0.b and huo /t<l 0. For T»W p /k B , 
the exponential prefactor in Eq. 3 will be proportional to 
T/uj ||. In this case, the exponent of the isotope effect 
on the prefactor is about -0.5, in sharp contrast to a large 
and positive exponent predicted by Eq. 3. 

In Fig. 3, we show the resistivity in the paramagnetic 
state for the 16 O and 18 O films of (a) LCMO; (b) NSMO. 
It is clear that the resistivity strongly depends on the 
oxygen isotope mass. The data can be well fitted by 
Eq. 3 with the parameters: A = 17.3/12.7 mftcm K°\ 
E p /k B = 844.3/997.2 K for the LCMO 16 0/ 18 film; 
A = 22.2/16.2 mficm K 5 , E p /k B = 925/1082 K for 
the NSMO 16 0/ 18 film. For LCMO, the parameter A 
decreases by 35(5)%, and E p increases by 13.2(5) meV. 
For NSMO, A decreases by 37(8)%, and E p increases by 
13.5(9) meV. The observed large oxygen isotope effect 
on the parameter A is consistent with Eq. 3. It is worth 
noting that since the activation energy is large, one can 
also fit the data by p(T) = AT 13 exp(E p /k B T) with (3 
varying from -2 to 2. 

We can use the values of the parameter A for the 16 O 
films to calculate the polaron bandwidth W p according 
to the relation: A = (ah/e 2 ^k^)(1.05W p ) 1 - 5 /hw . By 
taking hu> — 74 meV, typical for oxides |30|, we yield 
W p = 49 meV for the LCMO 16 O film, and 58 meV for 
the NSMO 16 film. The polaron bandwidth is greatly 
reduced compared with the bare bandwith W 0l which 
can be estimated to about 4.9 eV from the relation: 
W = 6h 2 /a 2 rrib and mi, = 0.61m e p8[ . Then we can 
determine g 2 to be 4.61/4.44 for the LCMO/NSMO 16 
film using g 2 = In (W /W p ). If Tioj decreases by 5.7% 
upon replacing 16 O by 18 O, then our calculation shows 



that the parameter A decreases by 35% for LCMO, and 
by 33% for the NSMO, in good agreement with the mea- 
sured values: 35(5)% for LCMO and 37(8)% for NSMO. 

It is important to clarify whether the small polarons 
are bound into localized pairs (bipolarons) or just to the 
impurity centers since this can place an essential experi- 
mental constraint on various CMR theories. If the small 
polarons are bound to the impurity centers, there will 
be no isotope effect on A [|o|. In this case, the iso- 
tope effect on E p only arises from the isotope shift of 
E a , as seen from the relation: E p = E a + A/2. Since 
E a = (r)E p /2)f(T) - W /12, only the quantity f(T) = 
[ta.nh(huj /4:k B T)]/(huj /4:k B T) may depend on the iso- 
tope mass if the temperature is not so high compared 
with tku /k B . If we take hu = 74 meV, we have f(T) = 
0.855 and 5f(T) = 0.0134 for T = 300 K, where S means 
a change due to the oxygen isotope substitution. Us- 
ing the relation: {-qE p /2)f{T) = W /12 + E p - A/2, and 
A/2 = E s = 15 meV for La .75Cao.25Mn0 3 @, we yield 
-qEp/2 = 0.54 eV. Then SE p = 5E a = (r]E p /2)Sf(T) = 
7.2 meV, which is about half the value observed. 
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FIG. 3. The resistivity in the paramagnetic state for 
the 16 and 18 films of (a) Lao.rsCao.asMnOs; (b) 
Ndo.7Sro.3MnC>3. The solid lines are the fitted curves by 
Eq. 3. Note that the fit is not so perfect for the NSMO 18 
film since the data are not very smooth around T = \ATc- 

On the other hand, one can quantitatively explain the 
isotope dependence of E p if small polarons form local- 
ized bipolarons. In this scenario, SA — —5W P p3|, so 
SA = 0.057g 2 W / p . From the deduced values for g 2 and 

12.9 meV for LCMO, 



W p above, we calculate that 5 A 



and 14.7 meV for NSMO. Then 6E P = SE a + 5 A/2 



3 



13.6 meV for LCMO, and 14.5 meV for NSMO. The cal- 
culated values are in excellent agreement with the ob- 
served values: 13.2(5) meV for LCMO and 13.5(9) meV 
for NSMO. In addition, if the small polarons are bound 
into small bipolarons, there will be a maximum in optical 
conductivity at the energy i/& = 2g 2 Twj D + W p + A [ p9| . 
Substituting g 2 = 4.61, hu = 74 meV, W p = 49 meV 
and A = 30 meV for the LCMO 16 O film, we obtain v h = 
0.76 eV, in excellent agreement with experimental results 
(i.e., v b ~ 0.70 eV for x = 0.3, and 0.80 eV for x = 0.1) 
& 

In summary, the oxygen isotope effects observed in 
high-quality epitaxial thin films of Lao.7sCao.25Mn03 
and Ndo.7Si'o.3Mn03 can be quantitatively explained by 
a scenario Q where the small polarons form localized 
bound pairs (bipolarons) in the paramagnetic state while 
the carriers in the ferromagnetic states are small polarons 
due to the competing exchange interaction of polaronic 
carriers with localized spins. The electron-phonon cou- 
pling strength g 2 in the paramagnetic state is about twice 
of that in the low temperature ferromagnetic state. This 
is possibly because a strong screening effect by mobile 
charge carriers reduces the electron-phonon interaction 
in the ferromagnetic state. The present results provide 
quantitative constraints on the basic physics of mangan- 
ites. 
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